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A NOTE ON THE LOCAL REGULARITY OF DISTRIBUTIONAL 
SOLUTIONS AND SUBSOLUTIONS OF SEMILINEAR ELLIPTIC 

SYSTEMS 

RAINER MANDEL 


Abstract. In this note we prove local regularity results for distributional solutions and 
subsolutions of semilinear elliptic systems such as 

L'^Uk = fkix,ui,...,UN) in R” (fc = I,...,A^) 

where Ai,..., Aat are of divergence-form and n > 2m. We show that distributional subsolu¬ 
tions are locally bounded from above if \fk{x,z)\ < C(1 + |z|^) for 1 < p < , k = 1,... ,N. 

Furthermore, regularity properties of subsolutions and improved versions for bounded sub¬ 
solutions are given. Even for /i = ... = Jn = 0 our results are new. 


1. Introduction 

The starting point of this paper is the following question: Being given a semilinear system 
of elliptic partial equations of the form 

(1) L^Uk = fk{x,Ui,...,UN) in (fc = l,...,A) 

on some open set hi c R"- what is the maximal regularity of an arbitrary distributional subsolu¬ 
tion or solution of ([T])? More specifically we are interested in assumptions on divergence-form 
operators Li,..., Lat of second order and nonlinearities fi,. .., which ensure that distribu¬ 
tional solutions or subsolutions (see ([3j),(jl]) for a definition) are locally bounded or bounded 
from above, respectively. 

The study of unbounded weak and distributional solutions was initiated about 50 years 
ago and it gave rise to several interesting methods and results. Let us try to give a short 
overview of the subject with a focus on second order equations (i.e. m = 1, A = 1) such as 

(2) - dw{A\/u) = f(x,u) in n 

where hi is a bounded domain in R*^, n > 3. We start with unbounded solutions that owe 
their existence to the roughness of the matrix function A which we will always assume to 
be positive definite and bounded on hi. In a fundamental paper Serrin mi provided explicit 
examples for matrix functions A and unbounded solutions of ([2]) for / = 0 such that the 
solutions he in Sobolev spaces with s < 2. Given the fact that local regularity 

results of De Giorgi-Nash-Moser type (see chapter 8 in 0) are valid for weak solutions lying 
in we see that different a priori assumptions on the regularity of the solution may lead 

to different kinds of solutions. In the works of Brezis [1] and Jin, Maz’ya, van Schaftingen [TOj 
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this issue was analyzed further. For instance, they proved that if the entries of A are only 
assumed to be continuous then solutions lying in do in general not possess better 

regularity properties whereas solutions in p > 1 he in for all q £ [l,oo). 

In particular these solutions are always bounded (in contrast to their gradients, see HU]). 
Under the assumption of Holder-continuous coefficients the latter result had previously been 
obtained by Hager and Ross [9]. For weak solutions in Brezis proved the 

regularity for matrix functions A with Dini-continuous entries, see Theorem 2 in |1]. All 
these results concern weak solutions whereas Bao and Zhang w studied regularity results 
for distributional solutions with Lipschitz continuous A. They showed that in this case 
distributional solutions of -div(AVM) = / have the ’’typical” regularity properties of elliptic 
problems and no pathological solutions like the ones mentioned above can exist. 


Another way of producing unbounded solutions of ([2]) is to consider nonlinearities f(x,z) 
that grow sufficiently fast to infinity as z tends to inhnity. This can be illustrated via the 
model equation - An = on a domain H c For p > ^ unbounded distributional solutions 

are given by the explicit formula u{x) := Cn,p\x - for some Cn^p > 0,a;o e H. More 

sophisticated unbounded distributional solutions of this equation in suitable domains H c 
with appropriate boundary conditions are due to Pacard [IU1|TU] and Mazzeo,Pacard [TU] for 
exponents p > Finally let us mention that similar constructions were performed in [T2| 
in the context of nonlinear Schrodinger equations on R” for p slightly larger than 

Let us now describe in which way our results contribute to the issue. We deal with 
subsolutions (Theorem [T]) and solutions (Theorem [2]) of the 2m-th order semilinear elliptic 
system ([T]) on an open set H c R".. We restrict our attention to the case n > 2m which, 
from the point of view of regularity theory, is more interesting. We hnd regularity properties 
of subsolutions which will be shown to be optimal in a general setting. A new feature of 
our approach is that these results can however be improved once we add some integrability 
assumption on the negative parts of the subsolutions. Furthermore, even in the easiest 
case of linear problems of second order equations (m = 1,N = 1,/ = 0) our results are 
new since the involved linear divergence-form differential operators Li,...,L]\f may have 
Lipschitz continuous but also less regular coefficient functions, see assumption (AI)q, below. 
In particular we can treat more general situations than in [21|U] where distributional solutions 
of -div(AV'n) = / were investigated under the assumption that A is locally Lipschitz. Our 
assumptions on the coefficient functions will be shown to be sharp in the sense that for slightly 
less regular coefficients our regularity results cannot hold any more in view of the pathological 
solutions found by Jin, Maz’ya, van Schaftingen [10] and Serrin HZI. In Remark [1] (c) this 
aspect will be explained in detail. 


Before coming to the statement of our main result let us provide the definitions of distri¬ 
butional solutions and subsolutions of O- To this end we introduce a class of differential 
operators which is suitable for the definition of distributional solutions lying in L[^^(r2;R^). 
We will always assume that the following hypothesis is satisfied: 
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{Al)a The differential operators Li,..., Ljq are given by Lkcf) := - div(AfcV(/>) for cf) e 

where the matrix functions Ai,... ,An are symmetric and positive definite with entries 
in for some a e [1, :;^)- 

The assumption 1 < a < ^ makes sure that Ai,... ,A]^ are Holder-continuous by Sobolev’s 
embedding theorem so that several classical results in the theory of elliptic partial differential 
equations [7] can be applied in the sequel. In Remark [T](c) we will show that this assumption 
is not only helpful from a technical point of view but also essential for our Theorem to be true. 
We say that u e is a distributional solution of ([I]) if we have f(-,u) € 

as well as 

(3) f UkL]^(j)= f for all 0 6 (k = l,...,N). 

Jn Jn 

Accordingly a distributional subsolution is supposed to satisfy 

(4) [ UkL^(j)< [ M;u)cP forall06C'r(H),0>O {k = l,...,N). 

Jn Jn 

Next let us state the main result of this paper. 

Theorem 1. Let H c M”' be open, n > 2m and assume u e solves 

L'^Uk < C{g + + ... + in Vt (k = 1,..., N) 

in the distributional sense with and where satisfy (Al)^ and 

g e L'[^J^VL),r > Furthermore assume 1 <p < or p> j ^ some 

q > ■ Then we have e and u e R-^) for all t e [1,^). In 

addition, the following implications hold true: 

(i) Ifn>2m and e g > then u e iy“’*(H;R^), I <i < ■ 

(ii) Ifn = 2m and u. e then u € (H;R^). 

Here we used the notation tt* := ... ,u^) to denote the vector of the positive/negative 

parts of the component functions of u. The proof of Theorem [T] is based on a representation 
formula for subsolutions and a well-known bootstrap procedure that seems to go back to 
Stampacchia in the case m = 1. In the following remark we discuss extensions of Theorem [T] 
and why it can not be essentially improved. 

Remark 1. 

(a) The iteration scheme from the proof may be slightly modified to prove local boundedness 
results for problems of the kind 

Lffuk<C{gi +g 2 uY‘+ ... +g 2 u]^^) inVt {k = l,...,N) 

where gi € ^ and g 2 e L[/^(H),r 2 > 1. In this situation the critical 

exponent for local boundedness changes to (1 - T), Also the differential 

operators L™ = LkO .. .o Lk may be replaced by compositions of m different divergence- 
form operators as in (Al)^ without changing the result. 
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(b) The restrictions on p, q from the theorem are optimal in view of several results for the 

model equation -Au = u^. Unbounded distributional solutions in the case p = were 
found by Pacard m and in the case p > the function u from the introduction may 
be taken. Note that u e only for 1 < q < . 

(c) The condition a € [1, ;^) from assumption (Ai)^, is sharp in the sense that Theorem\^ 
is not true for a e [^,oo). Indeed, in m Proposition 1.2 the authors constructed 
a cotinuous matrix function A and a weak (and hence distributional) solution u e 

c of the equation div(v4Vw) = 0 such that u i for any 

other p> 1. Having a look at the proof of this result (and in particular Lemma 2.1 and 
equation (8) in fWf) one realizes that the coefficient functions of A lie in 
Hence, this example shows that our theorem does not hold in the case a = For 
Oi £ Theorem can not hold either since Serrin’s pathological solutions 

from mi provide counterexamples. In the case a e the nonlinearity allows 

us to take the solutions from (b) as counterexamples. It is unclear to the author, 
however, whether the linear problem div( 74 VM) = 0 admits unbounded distributional 
solutions for such a. Finally, we think that the case a = oo has not been treated yet. 

(d) Also the results from the parts are (i),(ii) are close to optimal. First of all one should 

notice that regularity results in Sobolev spaces of order 2m or higher can in general not 
hold since the fundamental solution of (-A)™' does not have such weak derivatives. 
Moreover, the functions us(x) := a(xf + ... + for small 5 > 0 and appropriate 

choice of a € define bounded poly subharmonic functions which do not lie 

(M*^) so that our result for q = oo may be considered as sharp. 
More generally, for suitable a e the functions x + ... + x\yF define 

poly subharmonic functions in provided 6 > max{ 2 m - 

These functions do not lie in certain q the exponent 2 m-i-s 

can be very close to i_|_( 2 m-i)q - precisely, if we could define these functions for 

all k ■■= 6 (2m, n) (which, in general, is not a natural number) then we would 

obtain the optimality of the exponent. Unfortunately, we have to leave open whether 
non-formal examples exist or not. 

Theorem [T] admits a refined version for distributional solutions of ([T]) which we will formu¬ 
late below for the sake of completeness. It generalizes the results of Bao and Zhang to 
nonlinear higher order problems and complements the existence and regularity results of Jin, 
Maz’ya, van Schaftingen [10] and Brezis [1] in the sense of Remark [T] (c). The assumptions 
on the right hand side / are the following: 

(A2) f : Q X R^ is a Caratheodory function satisfying 

\f(x,z)\<C(g(x) + \z\^) on n X where (7 e L[^^(r 2 ), r e ( 1 , oo). 

Under the assumptions (Al)a,(A2) we show that in the case 1 < p < ^very distribu¬ 
tional solution u € Lf^yn-jR^) n is a bounded weak solution so that classical 

elliptic regularity results as in [7] are applicable. We will omit the proof since it results from 
discussing special cases in the proof of Theorem |2l see Remark |2] (b). 
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Theorem 2. Let c M" be open, n > 2m and let u € n solve 

(5) L^Uk = fki,x,Ui,...,UN) inn {k = l,...,N) 

in the distributional sense where Lk,fk satisfy (Al)a,(A2). Moreover assume l<p< 

or p> for some q > Then we have u € W^ff’^(n]R^). Moreover, 

A € C^(n;R^x^),f € C°°(n X R^-R^) implies u e 

As a consequence the assumptions of the theorem guarantee that unbounded distributional 
solutions of ([5]) can only exist for n > 2m, p > these solutions have to be 

searched in subspaces of L‘l^^(n) with p < q < while distributional solutions lying in 

higher Lebesgue spaces are automatically locally bounded. 

In Section [2] we provide an auxiliary Lemma needed for the proof of the parts (i),(ii) of 
Theorem [U In Section [3] we prove Theorem [T] using some results concerning Green’s functions 
which we provide in the Appendix. 

2. An auxiliary lemma 

For a given ball B c R'^, 7 < 0 and a Radon measure p on B with p{B) < 00 we introduce 
the measurable functions : R ^ [0, cx)] as follows: 


( 6 ) 


JB Jb 


By definition, a Radon measure on R is a Borel-regular measure on R which is locally hnite. 
In the proof of Theorem [1] we will need the following result. 

Lemma 1. Let n,k e N,j,p,q € M satisfy 0 > j > k - n,q > p > 1 and let fj,,B be given as 
above. Then the following implications hold true: 

(ti 'y^Q 

(i) If (f)^ e Li (B) then cf^-k e L^ (B) provided l<p< - - -—, 

n + 'y + {q-l)k 

IL 

(a) If 4> ^ Llf^B) then cf-k e Lf^^B) provided n > 2k. 

Proof. In order to prove part (i) we set r := so that the following is true: 


— > r > 1, 


kr 


7- 


> -n. 


p r - 1 

Using these inequalities and Holder’s inequality we obtain for every B B 

J_\(j)^.k{x)\P dx = ~ y\^ \x-y\^^^'" dfi{y)J dx 


P(r-l) 


Jjx - yl'f dfi{y)y Jjx - y\'^ y dn{y)'^ " 
J^(j)yx)r Jjx - y^^ dfi(yy " dx 
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= Ml^B)(f^fj^-yr^^My)dx 

qr—p 

= X( f^\^-y\^~^^d^)My)) " 

qr—p 

P. qr—p 

for some C > 0. This proves part (i). The proof of (ii) is similar. Using analogous estimates 
for 7 = 0,g = oo^p = we find that (p^k lies in Lp(B) provided ip € L°°(B) and 

sup I \x - y\~'^ log(l/|a: - y\y~'^ dx < oo. 
y^B Jb 

This condition is satisfied if and only if n > 2A: which proves part (ii) of the Lemma. □ 

Remark [1] (c) tells us that at least for some special choices of q this result can not be 
improved in an essential way. 


P(r-l) 


qr 

qr—p 


dx^ 


qr—p 

qr 


3. Proof of Theorem \T\ 


From now on we assume that the hypotheses of Theorem [T] are satisfied. In particular 
we will use the hypothesis {Ai)a where 1 < a < Without loss of generality, we may 
consider only the case = 1 which simplihes the notation in the following. First we discuss 
the regularity of subsolutions of linear equations. We use the following result which, for 
Lipschitz continuous matrix functions in the case m = 1, is due to Bao and Zhang. We only 
present the main changes with respect to the proof of Theorem 1.3 in [2] . 

Proposition 1 (Linear regularity I). Let h e L^^(r2) satisfy L'^h = 0 inQ in the distributional 
sense. Then h e for all t € [1, oo). In particular h e 


Proof. First we consider the case m = 1. So let /3 = > a and our hrst aim is to show 

h € To this end we proceed as in [2] and take an arbitrary function w e (7°°(hi), let 

R cc n be a ball and let (p be the uniquely determined function satisfying 

(7) -dw{AV(p)=w inR, 0 € Wo’^(R) n lU^’^(R), 

see Theorem 9.15 in [7]. By Theorem 9.19 in [7] the function cp is twice continuously dif¬ 
ferentiable with Holder-continuous second order derivatives and the differential equation is 
satished almost everywhere. Moreover, we have 


II0II 




(B) 


W •'^(B) 


<c'l 


w\ 
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for some C,C' >0 by Sobolev’s embedding theorem and Lemma 9.17 [7]. Therefore we may 
test the equation with rjcj) for some cut-off function r] e C^{B) and obtain 

0= f h{-dw{A\/{r](l)))) 

JB 

= I hi-div(A\/(f))r] - 2V(fFA'^rj - div(A\/r])(j)) 

JB 

L~(S)) 

- hwf] - C 

> f^hwr]- C" II h lU.(B) II71II ^ II ic II 

by definition of (3 for some (7, C', C" > 0 independent oiw. As in [2] the dual characterization 
of Lebesgue spaces gives h e and iterating this procedure as in [2] yields h e L\^J^B) 

for all f < oo. 

In order to prove the existence of weak derivatives we can not proceed as in [2] since 
there a difference quotient method is used which relies on the Lipschitz continuity of the 
matrix function. Instead we continue with the duality argument. To this end we consider 
Wo,Wi,... ,Wn ^ C°°{B) and choose (f) to be the uniquely determined function satisfying 


(8) -div(AV(/)) = + (f) ^Wl'‘^{B) nW‘^^‘^{B). 

i=i 

As above, 0 is twice continuously differentiable with Holder-continuous second order deriva¬ 
tives and the differential equation is satisfied almost everywhere. Instead of Lemma 9.17 
in [7] we use the estimate 

n 

II011 1 t , < C ||wd| t 

from Lemma 2.2 (i) [5] for t > > a. Proceeding as above we find (7, (7', C" > 0 such that 

0= f h(-div(Av(r](j)))) 

JB 

P n 

> / /i(u;o + E^i^i)7-C'II^IU7s)(IIII^IU~(s)ll|V0liy^.g) + ||A||^^i._^(^^||0||^) 

JB ^ 

p n n 

> h{wo + djw,)r] - C"\\h\\LtiB) ||A||^1._^ y; IIw,- 


In the second inequality we used Sobolev’s embedding theorem and a < From this 

estimate and the dual characterisation of see Theorem 3.9 in [1], we arrive at 

h 6 W^’^{B) for all t > As in [2] this induces h e W‘^^’^{B) for all t < oo. 
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In the case m > 2 instead of ([H]) (i.e. Lcj) = w in B) one solves L'^cp = ta in i? in some 
snbspace of which can be done by indnction over m nsing the same theorems 

as above. Similar estimates then allow to conclnde as in the case m = 1. □ 


With Proposition [T] at hand, we may now discnss the regnlarity properties of snbsolntions 
of linear problems. One main featnre of onr result is that integrability assumptions on the 
negative parts of subsolutions can be used to deduce slightly better regularity properties. 


Proposition 2 (Linear regularity II). Let u e satisfy L'^u < g in Q in the distribu¬ 

tional sense where g e L[^J^Vt),g > 0. Then the following implications hold true: 

(i) Ifr = l then e for all s e [1, 

(ii) Ifr^ [1, then e 
(in) Ifr>f)^ then e 

In each of these cases one has u e iP) for all t e [1,^). Moreover, assuming 

and r > 2 mnq+rff^{n- 2 m) following implications hold true: 


(^v) Ifq> ^,n> 2m then u e IP“’‘(0) for all t € [1, 
(v) If q = oo,n = 2m then u e pp^ 2 m-i, 2 m/( 2 m-i) 


Proof. It suffices to verify the above-mentioned regularity properties in an arbitrary com¬ 
pactly contained ball B cc O. For any given such ball let Gi be the Green’s function of 
the operator L on a slightly larger ball B' associated to homogeneous Dirichlet boundary 
conditions on dB'. Under our regularity assumption {AV)a the existence of Gi is guaranteed 
by Theorem 1.1 in [8]. We define the function Gm : S' x S' ^ [0, cx)] inductively via 

(9) Gk{x,y)= / G{x,z)Gk^i{z,y)dy {k = 2,...,m). 

JB' 

Then one has L'^Gmi'^y) - - ?/) in S in the distributional sense as well as 

for |q;| = 1,... ,2m-1 {x,y€B), 


( 10 ) 


G^{x,y) < 
Gm{x,y) > 


(G\x - y\^'^ ”■ 
\G\og{l/\x-y\) + G' 
(c\x - ypm-’T’ 
\clog(l/|x-i/|) -c' 


, if n > 2m, 
, if n = 2m 

if n > 2m 
if n = 2m 


(x,y 6 S), 
{x,ye S). 


Here, c, c', G, C' are positive numbers independent of x, i/ e S, a e Ng is a multiindex and 5 is 
the Dirac measure on M** centered at 0. In the Appendix we provide the references for these 
estimates. 

Step 1: A representation formula. The function vb '■= -u\b + jB^m{')y)g{.y)dy satisfies 
> 0 in S in the distributional sense. Using Theorem 2.17 in [IT] we obtain L'^vb = ^b 
in S where /i^ is a Radon measure. Defining Lb '■= -vb + JB^m{',y) dyB^y) we arrive at 

u\b = - / Gra{-,y)dyB{y)+ / G^i', y)giy) dy+hB where L™he = 0,/is € L^^^^(S). 

JB JB 

' -V-" '-V-" 

=-.ui =:U2 


( 11 ) 
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Moreover, + g)\B = + g)\B')\B = ^b' -I B implies ^b^B) < cx) since is a 

Radon measnre and thus locally finite. This property will be used when we apply Lemma [1] 
in Step 3. 

Step 2: Proof of (i), (ii), (in) - Integrability. By the Hardy-Littlewood-Sobolev Theorem 
(see m, Theorem 4.3) U 2 has the integrability properties which we claimed to hold for u 
in (i),(ii),(iii), respectively. By Proposition [1] we moreover have Hb e Hence, the 

inequality < M 2 + |^b| implies that lies in the same Lebesgue spaces as U 2 which is what 
we wanted to show. 


Step 3: Proof of (iv),(v) - Regularity. From (ITU]) we get Mi,M 2 ^ for all f e [1, :^). 

Hence, Proposition [T] tells us that u = -Ui + U 2 + Iib lies in the same spaces. Now let us 
additionally assume u~ € L^{B),g e for q,r as in the statement of the theorem. The 

assumption on r, the upper bounds for the derivatives of Gm from (ITOl) and the Hardy- 
Littlewood-Sobolev Theorem imply 


( 12 ) 


112 elR;,, -(B) 


w, 


2m-l, 


2mg 


l + (2m-l)g 


loc 


(B). 


Furthermore, the assumption u e implies U\ + |/is| e Li(B) via fllip and thus 

Ui € L(^^(B) by Proposition [H In the case n > 2m this implies (t> 2 m-n ^ B'(^^(B) where (f) 2 m-n 
was defined in (jUj). Indeed, using the lower bound for from fITUD we get 

ui(x)= / Gm(x,y)dgB(y) >c \x - dgB(y) = c(l) 2 m-n(x). 

JB J B 

From Lemma [D (i) we obtain e Lp(B) for p e [1, upper bound for 

the derivatives of Gm from (ITU]) imply 

(13) u, £ if P Hi. .JiT 11- )’ " ^ 2 ™. 

1 + {2m - l)q 

Similarly, in the case n = 2m the assumption u~ e L(fJ^VL) implies e L^fJ^B) so that 
Lemma [1] (ii) yields = 0i_2m ^ and thus 

(14) Wi^LPioc '^-"(B) ifg = oo, n = 2m. 

From (1T2|) ■ (lT3i) ■ (1141) and Proposition H] we obtain the assertion of the theorem. □ 


Remark 2. 

(a) Every (nonnegative) Radon measure g on B defines a distributional subsolution via 
the formula CB. 

(b) In the more special case L^u = g we have dgB(y) = g~(y)dy which improves the 
regularity properties of ui according to the integrability properties of g~. This is re¬ 
sponsible for the fact that solutions of L^u = g with g e L'(^^(Q),r > 1 are more regular 
than subsolutions. For elliptic equations with measure-valued right hand sides (such 
as subsolutions) this is in general not true. 


Proof of Theorem [I] 
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Applying Proposition |2] to g ■■= C{1 + we find that it is sufficient to prove «+ € 

In case 1 < p < the proposition yields e and thus g e for some 

qo e (p, )■ Using again Proposition |2] one inductively proves where 


Qk+l ^ 


mk 


np~2mqk: 

2^/c 

OO 


if 1 < ^ 

p 2m ’ 

•r 2* _ 

p 2m ’ 

if 

p 2m' 


The sequence (qk) increases due to qo > (because of go > p) and reaches +oo after 

hnitely many steps so that the assertion is proved in the case 1 < p < In the case 

n > 2m,p > assumption «+ € L‘^^^JVL),q > leads to the choice go := g and the 

same iteration as above gives the result. □ 


Appendix - On Green’s function 

Let us briefly recall why the estimates from ffTOjl hold if the assumption (Ai)^ is satished. 
First we note that due to this assumption the (2m - 2)-th order derivatives of the matrix 
functions Ai,... ,A]^ from the theorem are Holder-continuous by Sobolev’s embedding the¬ 
orem so that, generally speaking, L^-estimates and Schauder estimates for elliptic problems 
in divergence-form are applicable. In the following let A be one of these matrices and let Gm 
be the function dehned in ([ 2 ]). 

We hrst analyze the properties of Gi and we start with the exceptional case n = 2. The 
logarithmic bounds for Gi are proved in section 6 in [5]. Notice that for the lower bounds 
one uses that B is compactly contained in B'. For the upper bounds of the derivatives let us 
recall from [B] and estimate ( 2 ) in [S] that the following interior estimates hold for any weak 
solution w of div(AVrc) = 0 in B': 

(15) < Gd”^||Vw||L 2 .oo( 52 ^) (1 < fc < 2m - l,n = 2) 

Here, B 2 d is a ball of radius 2d contained in B' and L‘^’°°{B 2 d) denotes the Lorentz space 
(or weak-L^ space). We refer the interested reader to the paper of Dolzmann and Muller [5] 
where these spaces are used in the context of Green’s functions. From estimate (11) in [5] 
we get IIVGi(-,p) 12 , 2 , 00 ( 5 ,) < G so that the interior estimates flTSD applied to Gi{-,y) on the 
ball B 2 d{x) with d= | min{|x - p|, dist(H, <91?')} yield 

|9"Gi(x,p)| < G\x - {x,y e B, 1 < |a| < 2m - 1, n = 2), 

see also [5], Lemma 3 where the same reasoning was used. This proves the estimates for Gi 
from ffTnjl for n = 2 . 

In the case n > 3 the bounds for Gi{x,y) follow from [ 8 ], Theorem 1.1. The bounds for the 
derivatives of Gi{-,y) follow from interior estimates in the same manner as above. Instead 
of flT5|l . however, one uses 

||v''ic||L~(Sd) < Cd~'"\\w\\L<-‘(B2d.) 


(k < 2m - 1, n > 3) 
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This Schauder type estimate follows inductively from [7], Corollary 6.3, see also Lemma 3.1 
in |8] for the special case k = 1. Hence we obtain 

\d'^Gi{x,y)\ < Clx - {x,y e B, jal < 2m- 1, n > 3) 

so that the estimates for Gi in the case n > 3 are proved, too. 

Finally, using the estimates for Gi one can inductively derive the corresponding estimates 
for Gm via the formula ([9]). 
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